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The effective potentials of the rf-SQUID and three- Josephson junction loop with a penetrating 
external magnetic flux are studied. Using the periodic boundary condition for the phase evolution of 
the wave function of Cooper pairs, we obtain new periodic potentials with cusp barriers in contrast 
with the usual smooth double-well potential. The tunneling through the cusp barrier becomes 
dominant for a parameter regime where the self inductance of the superconducting loop and the 
Josephson coupling energy are large. Calculating the tunneling frequency we see that the cusp 
potential may offer a more efficient qubit. 
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I. INTRODUCTION 

The rf-SQUID has attracted great attention in connec- 
tion with the study of quantum two-level systems^. In 
particular it has offered prospects for macroscopic quan- 
tum phenomena, such as quantum tunneling and quan- 
tum coherence between two quantum states^. A fre- 
quently appearing structure is the well-known double- 
well potential with a smooth barrier in some parameter 
regime. Surprisingly, however, the inspection into the rf- 
SQUID with precise boundary conditions has not been 
achieved. 

A periodic boundary condition can be obtained by a 
gauge invariance and the phase evolution by the wave 
vector of Cooper pairs which is identical with the fluxoid 
quantization condition. Using this boundary condition, 
we can obtain the effective potential of an rf-SQUID and 
three- Josephson junction loop which consist of the en- 
ergy of moving Cooper pairs and the Josephson coupling 
energy. The external flux induces a persistent current 
of Cooper pairs in superconducting loop. When the self 
inductance of the superconducting loop and the Joseph- 
son coupling energy are large, the total flux, sum of the 
external flux and induced flux, become quantized, while 
the fluxoid quantization condition is always satisfied. 

Considering the energy of Cooper pairs, the effective 
potential of the system shows a new periodic structure 
consisting of cusp and smooth barriers. In one extreme, 
the tunneling through smooth-barrier potential is domi- 
nant, which is the usual ca.sel&2££&, while, in the other 
extreme, the cusp barrier potential is more appropriate 
in tunneling phenomena. Both may be used as a qubit 
in quantum computing 7,8 . We compare the tunneling 
rates for the usual smooth-barrier potential and the cusp- 
barrier potential and show that the latter may offer more 
effective qubit. 

II. RF-SQUID 

We first consider an rf-SQUID shown in Fig. (a). 
We derive the effective potential of the system showing 
double-well shapes with both smooth and cusp barriers. 



The current through the Josephson junction in Fig. ^( a ) 
is written as 

I=-I c sm(p+-+CV, (1) 
K 

where </> is the superconducting phase difference across 
the Josephson junction and I c , R, and C are critical cur- 
rent, resistance, and capacitance, respectively. Using the 
relation, / = (<f> t — Q cy ±)/ L s , and the Josephson voltage- 
phase relation, V — —(&o/2ir)(f), Eq. becomes 

$ t = $ ext -L 5 / c sm0-L s ^ (t+ C <n ■ ( 2 ) 

Here $t is the total flux, <j> ex t the external flux, <E>i n d = 
L S I the induced flux, L s the self inductance of the su- 
perconducting loop, and <f> = h/2e the unit flux quan- 
tum. When L s I c /Qq << 1, one can use the approxima- 
tion (j) w 27r<I> t /$o m Eq. J3J), which results in the usual 
smooth double-well potential^. 

The fluxoid quantization condition for a thin super- 
conducting loopS, 

- $t + {m c /q c ) <p v c ■ dl = n$ (3) 
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FIG. 1: (a) rf-SQUID and (b) three- Josephson junction loop 
with the Josephson coupling energy Ej's and the wave vector 
of Cooper pairs k. 
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with q c = 2e, m c = 2m e , and the average velocity 
of Cooper pairs v c , can be represented as the periodic 
boundary condition, 



k n l = 2mi + 27r/ t 



(4) 



determining the wave vector at a quantum state n, where 
k n is the wave vector of the Cooper pairs, I the circumfer- 
ence of the loop, and ft = $t/$o- Here, k n l is the phase 
accumulated along the circumference of the loopifi. For a 
rf-SQUID in Fig. 2] (a) the boundary condition with the 
phase difference across the Josephson junction 4> becomes 



k n I = 2irn +- 2irft 



(5) 



The current, on the other hand, is given by / = 
— (n c Aq c /m c )hk n , where n c is the Cooper pair density, 
A the cross section, and the induced flux <5> in d = L S I 
becomes 



$0 



2tt 



(0) 



where 7 = (q 2 n c A/lm c )L s . Combining the expression 
for the total flux, 



ft 



L.I 



with Eq. JSJ, the wave vector is written as 
(1 + 7 )fc„Z = 27rn + 2tt/ - <f>, 



(7) 



(8) 



where / = $ ex t/$o- 

Neglecting the dissipative term originating from the 
quasi-particle tunneling through the junction, Eq. l|T)l 
becomes 

- (n c Aq c /m c )hk n = -J c sin</> - C($ /2n)4>. (9) 
With the representation for k n in Eq. JSJ we get 



$0 
2tt 



C<f> 



(1+7)tt 



2tt 



sin 0,(10) 



where Sj = I c &q/2it is the Josephson coupling energy 
and eo = 2ir 2 Ti 2 n c A/m c l. Eq. (|10(l describes the motion 
of a particle with kinetic energy Eq = Q 2 /2C with Q = 
C($o/27r)0 in an effective potential U c g(4>), 



U eS ((f>) = E n + Ej(l - cos< 



where 



E n = 



eo 



1 + 7 



/- 



2vr 



(11) 



(12) 



The meaning of E n can be understood more clearly 
when we represent E n as E n = (1 + 7)£'kin using Eq. 
iJSJ, where i^kin = (h 2 k 2 /2m c )n c Al is the kinetic energy 
of Cooper pairs. Since the induced energy of the loop can 
be represented as E'ind = (1/2)L S I 2 = 7-Ekin, we can see 



that E n is the sum of the kinetic energy and the induced 
energ y 11 ; 12 , i.e., E n = E kin + E ind . If we introduce the 
kinetic inductance Lk = m c l/An c q 2 , the kinetic energy 
becomes -Ekin = (l/2)£fc/ 2 and the parameter eo in Eq. 
(|10fl is represented by eo = ^/2Lk- If we introduce the 
energy scale Ek = ®l/2Lk = e and El = ^l/2L s , the 
parameter 7 becomes 



L s 

7 = T k 



Ek 
E L 



and the energy E n in Eq. I|12l) can be rewritten as 

E n = —— [n + f - — 
1 + 7 \ 2n 



(13) 



(14) 



The similar expression was obtained for a superconduct- 
ing loop without Josephson junction^. 

The number of excess Cooper pair charges on Joseph- 
son junction iV c = Q/q c is conjugate to the phase differ- 
ence <p and the commutation relation [(f), N c ] = i gives the 
quantum phase fluctuations of the junction. The Hamil- 
tonian H = P 2 /2M + U c s((f>) describes the dynamics of 
the phase difference, where P = N c h = —ifid/d(f> and 
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FIG. 2: Lowest energy giving effective potential(solid line) in 
units of E.j for /=0.49 as a function of Josephson phase for (a) 
?7=10 and (b) 77=5. The dotted lines represent the Josephson 
coupling energy. 



M = ($ /27t) 2 C. The charging energy Eq causes quan- 
tum tunneling between two wells in the potential U e ff((j>) 
of the rf-SQUID and the following three- Josephson junc- 
tion loop. 

The double- well structures shown in Fig. [21 are given 
by taking the lowest energy for the quantum states n = 
— 1,0, and 1 for the effective potential U c g ((f)). The shape 
of the effective potential is determined by the parameter 
r\ defined by 

For large 7, the value of 77 can be determined by two 
parameters L s and Ej. We calculate the value of the 
parameter 7 to obtain 7 « 600 for the parameters of the 
rf-SQUID of Ref. [4]. At much low temperature T — 
O(10mK), we roughly estimate the Cooper pair density 
in a pure metal n c ss n/2 9 with the conduction electron 
density n. Even if we consider the effect reducing the 
Cooper pair density such as band structure of the mate- 
rials, finite temperature, and so on, we expect that the 
value of 7 is still large such that 7 > > 1 . In this case the 
parameter 77 is approximately represented as follows, 

Ej 2L s Ej L S I C 

where Ej = I c $ /2ir with critical current of the Joseph- 
son junction I c . We obtain that 77 8.5 for the rf- 
SQUID*. 

For the value of 77 = 10, we can see the double-well 
structure in Fig. (a) with a smooth barrier centered 
at </>/2tt = 0.5 for $ oxt /$ =0.49. The usual double- 
well structure with a smooth barrier corresponds to this 
one. The left (right) well corresponds to the diamag- 
netic (paramagnetic) current state. Thus the total flux, 
$t — ^cxt + ^ind, for the state at left (right) well ap- 
proaches 0(<i>o) as the self inductance becomes large. For 
a smaller 77, however, a new double- well potential emerges 
around the cusp as shown in Fig. [21 (b) . The high cusp 
barrier shrinks while the low smooth barrier grows up 
as 77 decreases. Thus we get two wells separated by a 
cusp barrier. The state of current at right (left) well is 
diamagnetic (paramagnetic). 

At the local minima of U e g (</>), the condition 
dU e s(<p) I dcf> = gives —(n c AH 2 /m c )k n + Ejsmcj) = 
and the current relation becomes 
n c Aq, 



- Ttkf 



I c sin0 = 0. 



(17) 



Actually Eq. Q is also a current relation, i.e., I — 
($t — $ ext )/L s . If one consider both <I> t and / as in- 
dependent dynamic variables while <E> ex t is a control vari- 
able and minimize the energy of the system with respect 
to the variable <&t, the current relation, Eq. Q, can be 
obtainedii*A£. Now the remaining dynamic variables are 
(j) and / or equivalently k n . The current relation, Eq. 
(|17|l , is obtained by minimizing the total energy with re- 
spect to the variable <p. This current relation warrants 
the validity of our effective potential. 




FIG. 3: The effective potential of three- Josephson junction 
loop in (<j)a,4>b) plane, where Ej3/Ej=0.8, /=0.495, and 
77=1000. We set c /27T=//3. 



III. THREE-JOSEPHSON JUNCTION LOOP 

Recently the superconducting loop with three- 
Josephson junctions of Fig. d(b) has been proposed as an 
efficient qubil£. The periodic boundary condition similar 
to Eq. ||8J is written as 

(1 + 7)W = 27m + 2tt/ -0i -02 (18) 

The effective potential U e g((j)i,4>2,4'3) f° r this qubit is 
given by 

U eS (<t>i,<f>2,<f>3) = E n + Ej(l - cos^i) +Ej(l - cos0 2 ) 
+ B J3 (l-cos03), (19) 



where 



,-. E K I p "I-"-! - ' >:\ 

E n = I n + J 



1 + 7 



2vr 



(20) 



Here we introduced the rotated coordinates ((j) a ,(j)b,4>c) 
such that a ee (l/2)(0 1 +0 2 ),0 b ee (-l/2)(0!-0 2 ), and 
<j> c ee (l/3)(0i + (f>2 + 4>3). Then the effective potential of 
this qubit is written as, 

E ( 3 
^cff = — ^ (n + f - —(, 
1 + 7 \ 2n 

+2Ej{\ - cos (f> a cos 4> b ) + E J3 [1 - cos(30 c - 2<f> a )} . 

(21) 

Fig. [3 shows the effective potential at 4> c /2tt = //3, 
where we can see that the value of 4>b is strongly local- 
ized around <pb = nir since the barrier between Q and 
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FIG. 4: Contour plots for the effective potential of three- 
Josephson junction qubit with (a)?7=50 and (b)r;=2 when 
/=0.495 and (j>b=0. A, B, C, and D denote the ground states 
at the local potential minima. The lines <f> c w in (a) and 
(b) represent the cusp barrier. The potential in (c) is drawn 
along the line connecting the points A and C for rj = 50, while 
the potential in (d) is for the points A and D for r\ — 2. 



R is much higher than that between P and Q. This 
means that the phase differences at the potential min- 
ima in identical Josephson junctions are 2ivn. For 4>b = 
the effective potential is shown in Fig. 0] The barriers 
between A and C as well as between B and D are smooth, 
as can be seen in Fig. 0(c). The two current states at A 
and C correspond to the qubit states of Ref. [5] , where 77 
is so large that the cusp barrier is high enough for tunnel- 
ing across the barrier to be neglected. The cusp barrier 
between A and D shrinks while the smooth barrier grows 
as 77 decreases. In Fig. 01(d), thus, we can see the cusp 
barrier between A and D. The tunneling through this 
cusp barrier becomes dominant as 77 becomes small. 

The current in the loop is given by I = —I c wa.<f> a . 
Thus the current states at A and B (C and D) are in a 
paramagnetic (diamagnetic) current state. The values of 
4> a and 4>c at potential minima can be obtained through 
the current relations, dU c s/d<j>i — which give 



E K 



n + f- 



h + 4>2 + 4>3 



EjiShifa, (22) 



(1 +7)7T V" ' ' 2tt 
where Ej\ = Ej2 — Ej. In the rotated coordinate we 



sin < 



Asin(30 c -20 a ), (23) 



where A = Ej 3 /Ej and we set 4>b = 0. From Eqs. ltTH|) 
and (|23|l . we calculate / = (1 + -f)k n l/2ir. We calculate 
the values of the parameter 7 to obtain 7 w 200 and 
77 w 3300 for the three- Josephson junction qubit of Ref. 
[5]. For 7 >> 1, the induced flux in Eq. © can be 
represented such that ^md/'-'o = —lk n l/2ir w — /. 

In Fig. the induced flux is shown when / = 0.45 
and A = 0.8. The ground states in C and D induce a 
diamagnetic current while the ground states in A and B 
induce a paramagnetic current. In the inset of Fig. 
we can see that / is very small for the qubit of Ref. [5] 
so that the induced flux is negligible and the boundary 
condition, Eq. Q18f) . becomes 27m + 27r/ = (</>i + $2 + 
(^3) which is used in Ref. [5]. On the contrary, when 
the value of 77 is much smaller, the induced flux become 
larger. The value of total flux approaches ($0) f° r 
diamagnetic (paramagnetic) current, which clearly shows 
the flux quantization condition. Thus we can see that the 
total flux becomes quantized as L s Ej increases, while 
the fluxoid quantization condition, Eq. Q, is always 
satisfied. For superconducting ring without Josephson 
junction, the total flux is quantized when L s is sufficiently 
largeii. 



IV. TUNNELING FREQUENCY 

The quantum states in a cusp potential may be used 
to make a more efficient qubit because of high tunnel- 
ing frequency across the barrier. This allows more single 
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FIG. 5: Induced diamagnetic (solid line) and paramagnetic 
(dashed line) fluxes as a function of r\ when /=0.45. Inset 
shows the induced diamagnetic flux for large r\. 
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FIG. 6: Tunneling amplitude for (a) cusp barrier and (b) 
smooth barrier in the effective potential of rf-SQUID, where 
we set a=64000 (8000) and £j=300K (75K) for the cusp 
(smooth) barrier. Figures (c) and (d) are for the three- 
Josephson junction qubit, where we set q=64000 (80) and 
Ej=300K (10K) for the cusp (smooth) barrier. 



qubit operations in a short dephasing time. The geomet- 
ric configuration of the cusp potential, such as its higher 
curvature and smaller width, causes higher tunneling fre- 
quency. We compare the tunneling rate calculated by the 
WKB approximation, the instanton method^, and the 
Fourier Grid Hamiltonian Method(FGHM)^ in Fig. © 
The instanton method calculates the tunneling rate t as 
follows, 



/ = h^uCnJ |^ exp 



(24) 



where S inst = M J^dr0 2 (r), M = ($ /2vr) 2 C = 

h 2 /8E C , uj — y/k/M is the frequency of small oscilla- 
tion in a well with k = d 2 U e fi{(j)) / d(j) 2 \ l p Tliili , and C is a 
factor which depends on the shape of the double- well po- 
tential. The instanton action is evaluated in such a way 



that Sinst = V2M JT^ ln \/U c s(4>)d(j), where -</> m in and 
^min represent two potential minima. The instanton sat- 
isfies the equation of motion (l/2)A/0 2 (r) = U e s[(j)(r)]. 
Integrating the equation of motion to obtain the instan- 
ton trajectory and comparing it with the asymptotic be- 
havior, Co is given by 



C = 4>n 



2Mu 



exp[w I], 



(25) 



where I is a constant obtained numerically. 

The values of Eq are chosen through the parameter 
a = Ej/Ec which is roughly proportional to 1/rj in 



experiments^^. In Fig. we see that the frequency 
of tunneling through the cusp barrier for low rj is much 
higher than that for the smooth barrier in both rf-SQUID 
and the three-Josephson junction qubit. The tunneling 
through the potential barriers in the effective potential, 
UcS, introduces the tunneling term in the single qubit 
Hamiltonian, 



H, 



qubit 



(26) 



where F(<J> oxt ) is the energy level separation between two 
current states and a z ^ are the Pauli matrices. The first 
(second) term in Eq. (|26|) causes the single qubit rota- 
tion, R Z {R X ). The operation time r op for the R x gate is 
inversely proportional to the tunneling frequency. Since 
the quality factor is Q$ = T^jr ori with the dephasing 
time T0, the present results show that the cusped barrier 
allows the quality factor that is an order of magnitude 
larger than occur with the smooth barrier. 

Lower r\ requires higher Ej as shown in Eq. I|16|l . 
Thus, a higher inductive flux in Fig. [SJcould be detrimen- 
tal to the phase coherence of the Cooper pair wave func- 
tion. However, recent experiments for a three-Josephson 
junction loop with a Cooper-pair island^ and a current- 
biased Josephson junction qubit 17 show that using a 
large Josephson junction the high quality factor can be 
achieved. A current-biased qubit of dc SQUID type with 
large Josephson junctions can be described by the present 
theory for low value of rj. 

The quantum phase slip may be relevant in a 
thin superconducting loop when the energy levels are 
degenerate^. The phase slip process causes the change of 
tunneling frequency and, thus, induces the decoherence. 
For the lower value of rj in Eq. 1)16(1 the loop inductance 
is larger, which means the loop is thicker. Therefore the 
qubit using the cusp barrier with the lower value of 77 as 
shown in Fig. [2 may have a longer dephasing time than 
that using the smooth one. 



V. SUMMARY 

Using the boundary condition for the phase evolution 
of the wave function of the Cooper pairs, we obtain a 
new double-well potential with cusp barriers for both rf- 
SQUID and three-Josephson junction loop. At the min- 
imum of the effective potential we show that the current 
relation is satisfied. This picture allow us to calculate the 
wave vector of Cooper pairs, fc„, as a function of r\ and, 
thus, obtain the induced flux to show that the flux quan- 
tization condition is satisfied when the value of rj becomes 
small. Moreover we calculate the tunneling frequency as 
a function of 77 which is the characteristic parameter re- 
sponsible for the shape of the potential. The tunneling 
frequency in the cusp potential is much higher than that 
in the smooth potential, thus, the cusp potential offers a 
more efficient qubit. 
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